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ABSTRACT 

In time-distance helioseismology, the time signals (Doppler shifts) at two 
points on the solar surface, separated by a fixed angular distance are cross- 
correlated, and this leads to a wave packet signal. Accurately measuring the 
travel times of these wave packets is crucial for inferring the sub-surface proper- 
ties in the Sun. The observed signal is quite noisy, and to improve the signal-to- 
noise ratio and make the cross-correlation more robust, the temporal oscillation 
signal is phase-speed filtered at the two points in order to select waves that travel 
a fixed horizontal distance. Hence a new formula to estimate the travel times is 
derived, in the presence of a phase speed filter, and it includes both the radial 
and horizontal component of the oscillation displacement signal. It generalizes 
the previously used Gabor wavelet that was derived without a phase speed filter 
and included only the radial component of the displacement. This is important 
since it will be consistent with the observed cross-correlation that is computed 
using a phase speed filter, and also it accounts for both the components of the 
displacement. The new formula depends on the location of the two points on 
the solar surface that are being cross correlated and accounts for the travel time 
shifts at different locations on the solar surface. 

Subject headings: Sun: oscillations; Sun: phase speed; Sun: spherical harmonics 



1. Introduction 



Time-distance helioseismology (iDuvall et al.l 119931 ) constructs wave packets by cross- 
correlating time signals between any two points separated by a fixed horizontal angular 
distance on the solar surface. It then measures their travel time bet ween the two points, 
by fit ting a Gabor wavelet to the observed temporal cross-correlatio n (IKosovichev fc Duvall 



1997 ). This travel time is then inverted to study different properties (IKosovichev 



199 



J 



Zhao 



20041 ) (e.g. sound speed perturbations, subsurface flows, meridional circulation) in the solar 
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interior that influence the wave packet, which cannot otherwise be studied by the global 
oscillations. The success in applying the technique depends on being able to design a phase 
speed filter of a certain width, that selects acoustic waves w ithin a certain rang e of horizontal 



phase speeds, and this improves the signal-to- noise ratio (iDuvall et al.l 119971 ) . Waves with 
the same horizontal phase speed travel the same horizontal distance on the solar surface, 
and therefore can be collectively used to probe the sub-surface features in the Sun, since 
they sample the same vertical depth inside the Sun. The filtering operation is carried out 
in the frequency domain, and it improves the signal-to-noise ratio by removing unwanted 
signals and waves that deviate a lot from the chosen phase speed, as these do not contribute 
to the cross-correlation, but instead can degrade it, and make the estimation of travel times 
inaccurate. The phase speed filter has a Gaussian shape centered at the desired phase speed, 
which is chosen to study waves traveling a particular horizontal distance. 

Previous approac hes to measuring the travel time between any two points was by fit- 



ting a Gabor wavelet (IKosovichev fc Duvalll 119971 ) to the measured cross correlation. This 



wavelet was derived by assuming that the amplitude of the solar oscillations has a Gaussian 
envelope in frequency of a certain width and peaked at a frequency where the power of 
the solar p-modes is concentrated, and moreover it considered only the radial component of 
the displacement. However, during data processing a phase speed filter is used. Also, the 
observed displacement on the solar surface has both radial and horizontal components. The 
horizontal component usually is ignored in the travel-time measurements. We show that 
it has a significant effect, particularly, for measurements far from the disk center and for 
moderate horizontal distances, contributing to systematic errors. 

To remove the shortcomings, we derive a new analytical cross-correlation wavelet that 
incorporates the phase speed filter, and also includes both the radial and horizontal com- 
ponents of the Doppler velocity oscillation signal. This wavelet retains the structure of the 
Gabor wavelet, and is a function of the filter parameters: central phase speed and width, and 
also on the amplitude, phase and group travel times that depend on the oscillation proper- 
ties and the dispersive nature of the solar medium. By including the horizontal component 
in the cross-correlation we also see a dependency on the location of the two points being 
cross-correlated, and moreover, the cross-correlation due to the horizontal component is a 
weighted sum of phase speed filtered Gabor wavelets, with the weights depending on the 
location of the cross-correlated points, the horizontal travel distance and the angular degree. 



Comparing this with the original Gabor wavelet formula (IKosovichev fc Duvalll 119971 ) 
we can estimate how the phase-speed filter and the horizontal component shift the measured 
travel times. For the line depth and intensity observations, that only have a scalar (radial) 
component, in the weakly dispersive limit which is true for wave packets that probe the deeper 
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layers in the Sun and travel large horizontal distances on the solar surface, corresponding to 
small values of angular degree both formulae are similar. Hence, using the old Gabor wavelet 
in this regime should not effect the travel time measurements due to the phase speed filtering. 
On the other hand, wave packets constructed by cross-correlating time signals separated by 
a small horizontal distance, are distorted by the dispersive medium in the outer layers of the 
Sun, and, hence, the new formula should be used to account for the travel time shifts due to 
the phase-speed filtering procedure. 



2. Cross-correlation of Phase Speed Filtered Signals 
2.1. Scalar intensity and line depth filtered signal 



In this section, we generalize the fitting formula of (iKosovichev fc Duvalll 119971 ) by 
including the phase-speed filtering. We first consider acoustic waves, that are observed by 
measuring either intensity fluctuations or line depth observations on the solar surface. They 
are represented as a sum of standing waves or normal modes at a point r = (r, 6, 0) in the 
solar interior and time t, and can be written as 

h{r, e, 0, t) = ^ anlmlpnlmir, 9, 0) exp {i{uJnlmt - Oinlm)} (1) 
n,l,m 

where, i = a/— 1, and each normal mode is specified by a 3-tuple {l,m,n) of integer param- 
eters, corresponding angular frequency oUnim, the mode amplitude anim, the phase anim and 
the spatial eigenfunction ipnim{f^i d, 4>) as a function of the radial variable r and the angular 
variables (0,(1)). The integer / denotes the degree and m the azimuthal order, —l<m<l, 
of the spherical harmonic 

Y^e, 0) = Q^P,"^(cos^^) exp(zm0), (2) 

which is a function of the co-latitude 6 and longitude 0. Here, is a Legendre function, and 
cim is a normalization constant. These describe the angular structure of the eigenf unctions. 
The third integer n of the 3-tuple {l,m,n) is called the radial order. 

For a spherically symmetric Sun the eigenfunctions ijjnir nir, 6, 0) can be s eparated into 



a radial function Xn«(^) and an angular component Y^^iO, (p) (jUnno et al.lll989l ). This repre- 



sentation is valid for example for the scalar intensity observations, 

'4Jnim{r,9,^) = Xni{r)Yr{9,(P) (3) 

and also all modes with the same n and I have the same eigenfrequency ujni, regardless of 
the value of m. In reality the Sun is not spherically symmetric, that causes this degeneracy 
in m to be broken. 
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In time- distance helioseismology we measure the travel time of wave packets by forming 
a temporal cross-correlation between the oscillation signals at two locations separated by an 
angular distance on the solar surface r = R, where R is the solar radius. To model this 
we represent the solar oscillations on the solar surface as a linear superposition of normal 
modes, that are band-limited in angular frequency u. 

H{R, 9, 0, Uj) = 271^ anlmXnl{R)Yl^{9, 0) exp{-ianlm)S{uJ - UJnlm) (4) 



where H{R,6,(j),uj) is the temporal Fourier transform (IBracewelll Il999l ) of equation (1). We 



ignore mode damping and assume spherical symmetry, he nce we can replace a;„;m by a;„/ and 



invoke the band-limited nature of the solar spectrum, (IKosovichev fc Duvalll 119971 : [Giles 



19991 ). by setting a„;mXnz(-R) = Gi{^), and we obtain from equation (jlj) the frequency band- 



limited signal Hg{R, 9, 0, u) = H{R, 9, 0, u)Gi{u!) 

Hg{R, 9, 0, a;) = 27r ^ Giiuj)Yri9, 0) expi-ia^iM^ - oj^i) (5) 

n,l,m 

where, the Gaussian frequency function Gi{uj) captures the band- limited nature of the am- 
plitude of the solar modes, and is given by 

G.H = &.exp(-^^^^) (6) 

The band-limited nature of the modes is controlled by the width 6uj, and a central frequency 
LJo, where the power of the modes is peaked in the u — I diagram. These two parameters 
change for different data sets, that probe different regions of the solar interior. 

A phase speed filter is applied to the signal to select modes from the u — I diagram 
for the purpose of constructing the cross-correlation wave packet. The phase speed filter is 
specified by a Gaussian centered around a phase speed Vph and a width 6Vph as parameters, 
and is given by 



(Vp - V, 



\2 



FM) = exp I - "^^^r 1 (7) 

Where, the phase speed Vp = j^, L = + 1) = k^R, is the horizontal wave number. 
The role of the phase speed filter is to select waves with a small range of phase speeds, 
the range is specified by the width 5Vph- AH these waves travel almost the same horizontal 
distance on the solar surface, and sample approximately the same vertical depth in the solar 
interior. Hence, it is crucial to select 5Vph appropriately so as to increase the signal-to-noise 
ratio in the cross-correlation function, and hence be able to better resolve the sub-surface 
structures in the Sun. 
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Due to the band-limited nature of the oscillation amplitudes, only values of L which 
are close to = ^ contribute to the sum in equation ([5]), and hence we Taylor expand L 
about the central frequency uJo, upto first order: 

L = L{uj) = L[ujo + (uj - UJo)] ~ L{i^o) + - ^o) (8) 

auj 

The equation ([8]) can be written in terms of the angular group velocity Ug = ^ and angular 
phase velocity Vp = evaluated at u = Uo, and using the fact L{uo) = y-, we have 

Likewise, the phase velocity Vp{L,u) can be expanded about the point {Lo,uJo) in the 
uj — I diagram to yield 



Fp{L,uj) ^ exp 



/ T/2 /^r _ JiL.V\ 



(10) 



where, 6f = iii^^^ and the filter width 6Vph is evaluated at (Lq, uJo), and is a constant. 

Phase speed filtering takes place in the frequency domain and consists of multiplying 
the filter function Fp{L,uj) with the band-limited Fourier transformed signal Hg{R,6,(f),uj) 
from equation 

Hf^{R, e, 0, uj) = 2nY, ^)Giiuj)Yrie, 0) exp(-m„;^)5(a; - lu^i) (11) 

n,l,m 

Equation (fTTj) can be inverse Fourier transformed to yield the phase speed filtered temporal 
signal hf^{R,9,(j),t) at a point R = {R,9,(l)) on the solar surface. This signal now consists 
of a superposition of modes that lie in a region of the u — I diagram, that is the intersection 
region of the frequency band-limited Gaussian envelope of the solar modes and the phase 
speed filter. 

hf^{R,t) = ^ Fp{L,u;ni)Gi{uni)Yi"'{9, (p) exp{-ianim) exp{iujnit) (12) 

n,l,m 

The temporal cross-correlation function il)f^{A,T) of the phase speed filtered signal 
hf^{R,9,(l),t) between two points A and B having coordinates Ri = {R,9i,(f)i) and R2 = 
{R, 92, 02) respectively on the solar surface with a fixed angular separation distance A, as a 
function of the time lag r is 

T 

^/,(A,r) = hj^iR„t)hf^iR2,t + T)dt (13) 
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where, hf^ is the complex conjugate of the phase speed filtered signal hf^, and T is the length 
of the time series. 

Substituting the expression for hf^{R,t) from equation ([T^ into equation ([T5|) one gets 
after applying the orthonormality of exp{—iujnit) in the temporal integral in equation (fT3l) 



tpf^{A, t) = Y^ F^{L, uJni)G^{uJni) exp{-iuJniT) ^""(^1' <Pi)Yr'{02, ^2) exp {-i{anim - anlm')} 



n,l 



m m' 



(14) 

where, Yl^ (6*2, ^2) is the complex conjugate of Yj^ {62, 4>2)- Since the phases are random, we 
assume that the term exp {—i{anim — anZm')} the double sum for m and m' is zero except 
when m = m'. In this case equation (14) becomes 



^;^(A,r) = ^F^2^L,.;„OG?KOexp(^^„,r)5^1^-(0i,</>i)>^r(^2,</'2) 



(15) 



n,l 



We can simpli fy equation (15) by applying the addition theorem of spherical harmonics 
(IjacksonlEgggj ) 



2tx 



(16) 



wher e, Pi is the Legendre polynomial of order I, L = + 1) ~ /+| for large Z ( Christensen-Dalsgaard 



20031 ). and A is the angular distance on the solar surface between the two points A (^1, 0i) 



and B (^2,02), and is (see Appendix A equations (H9l) and (150!) . Appendix D, Figure [6] and 
Appendix E, Figure [HI) 



cos A = cos 61 cos 62 + sin 61 sin 62 cos(0i — 02 ) 



(17) 



We can approximate for large /, small A, such that LA is large (I Jackson! Il999l ) 

■A~ 



Pi (cos A) ^ Jo 



{21 + 1) sin 



2 /r A 

cos LA 

ttLA V 4 



'1^ 



where, Jq is the Bessel function of the first kind. With this approximation from equation 
f fTSj) . the addition theorem in equation f|T6|) becomes 



5^rr(^i,0i)>l"^(^2,02) 



21/^ /rA ^ 

— W — COS LA 

/;^27r V 2 V 4 



(19) 



Using equation ( |T9i) . the equation ( |T5i) for the phase speed filtered cross-correlation becomes 
after taking the real part 



i^fpiA,T) =Y-^^^F^iL,UJnl)G^{Unl)cOs{UnlT) COS (^LA - ^ 



(20) 
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The double sum in equation fl20|) can be reduced to a convenient sum of integrals if we 
regroup the mo des so that the outer sum i s over the p hase speed Vp = and the inner 



sum is over Uni (IKosovichev &: Duvalllll997l : iGiles Ill999l ). The expression for ipf^{A, r) is an 
even function of r, so we get an identical expression if we replace r by — r (negative time 
lag). The product of cosines in equation (!20|) can be transformed by using the trigonometric 
identity, and this results in 



1 1 L 



COS ( UJnlT — LA 



+ COS ( CUnlT + LA — 

(21) 



In equation fl2T|) we let f+{uJniT) = cos [uJniT — LA + f ) corresponding to the positive time 
lag r and /-(cUn^r) = cos (cj^/t + LA — |) corresponding to the negative time lag — r. Since 
cosine is an even function, it is seen that f-iunir) = f+{—uJniT), hence we can drop f-iunir) 
from equation fl^ . while extending the sum to negative values of Uni- We now substitute 
the expressions for the Gaussian envelope frequency function and the phase speed filter from 
equations ([6]) and (fTOj) respectively into equation (!2T!) and since ujni is a dummy variable in 
the inner sum, it is replaced by uj for notational convenience. Moreover, the inner sum over 
uOni can be replaced by an integral over uj, hence the discrete inner sum in equation 0211) 
becomes 



oo 



\ 



exp 



COS {ujT — LA + 



(22) 

The integral in equation (l22l) is negligible for large frequencies hence there is very little error 
made in extending the frequency limit to oo and — oo. Here, we choose the coefficient hi for 
the Gaussian envelope function, such that 



(23) 



In order to evaluate the integral in equation fl22|) we Taylor expand L, and using the 
linear dispersion relation from equation ([9]) we get 



L 



UJ 
Vr 



ph 



UJ 



{ujT - LA) ^ UJ ( 



1 

T — Tn 



UJoiTp 



(24) 
(25) 



where, the group travel time is defined as = ^ and phase travel time Tp = ^ 
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Substituting equations flM|) and fl25l) into the exponent of the phase speed filter and 
the argument of the cosine term respectively in equation fl22|) . and carrying out the integral 
( IGradshteyn and Ryzhilall994j ). which is bounded since the exponentials decay at both the 
limits of integration, with some algebra yields 



V'/p(A,r, l/p) = Af^{6u;,6f,T,Tg,Tp) cos 



^ai'^g '^p) _|_ ^2 

Tph 



Where, Rg 



Tph 



and R„ 



Tph 



(26) 

are dimensionless quantities representing the relative 



deviation of the group and phase travel times respectively from the filter phase travel time. 



Tph 



T#- and 



ph 



A 



fp 



5uj e 



exp 



2^2 



(r - Tgf + 



IQujlRl 



(27) 



The argument of the cosine term can be written as 



Tph 



Tp) 2 \ ^o[^ - 'Tg j 

' Rl + e' 



(28) 



where, the shifted frequency, Uf^ = uJo{l — Rgps) and the shifted phase travel time due to the 
phase speed filter is 

^'"^ iRg-Rp)T,, (29) 



1-R 



gpe 



The dimensionless parameter Rgpe = -j^i^ is zero in the absence of phase speed filtering, 
which occurs when the filter width tends to infinity or in the non-dispersive limit when 



Rg = Rp = 0. Also from equation fl29|) we see a shift in the phase travel times due to the 
phase speed filtering, which tends to zero when Rgp^ tends to zero. Also if we choose Vph to 
be either Vp or Ug, then the effect of phase speed filtering can be removed since Rg — Rp — 0. 

Summing equation (l26ll over phase velocities we get the final cross-correlation. 

1 



V'/p(A,r) 



'TT 



A 



^f^{A,T,Vp) 



(30) 



In the non-dispersive limit when R„ = Rp = 0, that is t„ = r„ = Tph, or in the limit of no 



phase speed filtering, large e, equations 
formula 



exp 



phj 

and (1271) reduce to the original Gabor wavelet 



ir-TgY 



cos 



i^oiT - Tp) + 



(31) 
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In previous papers (iKosovichev fc Duvalll 119971 : iGiles Ill999l ). the factor of j has been ab- 



sorbed in the phase travel time. By defining the phase travel time r 



Tr, 



4a;o' 



the 



Gabor wavelet in equation flHTl) reduces to the familiar form with the sinusoidal part being 
cos [^o(r-rp)]. 

We conclude that while the phase-speed filtering procedure does not change the func- 
tional form of the basic time-distance helioseismology fitting formula, it systematically shifts 
the travel times, if the filter parameter, Vph, is different from the actual phase or group 
speeds for a given distance. 



2.2. Effect of Line-of-Sight Doppler Velocity Measurements 



We now consider the Doppler shift observations which are commonly used in helioseis- 
mology. The Doppler velocity depends on both the radial and horizontal components of the 
displacement. For simplicity, we take the axis of the spherical harmonics to be in the plane 
of the sky, orthogonal to the line of sight. 

The displacement vector d(f,t) at location r = (r,6, (t)) and time t is written for the 
spherically symmetric Sun as (jChristensen-Dalsgaardl 120021 ) 

'dY/^^ 1 dYf^ ^ 



n,Lm 



de 



sin 9 



exp {i{ujnit - anlm)} 
(32) 

where, ^^;(?"), ^nii^) are the radial and horizontal components of the displacement eigenfunc- 
tions respectively, and e^, eg and are unit vectors in r, 6 and (p directions respectively (see 
Appendix A, Figure E]). 

Thus, the observed Doppler signal on the solar surface is the projection of the displace- 
ment vector onto the line-of-sight direction. Without loss of generality we take the x-axis as 
the line-of-sight direction, hence the projection is just d{R, t) ■ e^, where e^: is the unit vector 
along the x-axis and is given by (see Appendix A, equation fH^ ) 



= ir sin 6 cos (t> + ig cos 6 cos 6 — 



, sm( 



(33) 



The Doppler signal on the solar surface is obtained by taking the dot product between d{R, t) 
and the line of sight direction e^, and is therefore d^{R,t) = d{R,t) ■ = '^nim^nimi^y'^)- 
For a mode (/, m, n) it is 



sin 6* cos + Pni{R) 



dY/^i 



09 



■ cos 6 cos (h — 



sin 6 dYi^' 



sin 9 



(34)' 
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where, 

^ " " ^ ^ ^ " ^ (35) 

The ratio of horizontal to radial component (3ni{R) shows dependence in frequency, and 
since the solar modes are band-limited in frequency, and peaked around cUo, we replace the 
frequency Uni by a constant uio = 27rz/o. Here, g is the acceleration due to gravity at the solar 
surface, a is the dimensionless frequency and be is a factor that depends on the boundary 
condition used at the solar surface r = R. For the solar 5-min oscillations = ^ ^ 1000, 
and for the boundary condition that the Lagrangian pressure perturbation vanishes at the 
solar surface: at r = i?, Sp = 0, leads to 6c = 1, and this shows that at the solar surface 
the radial component dominates. For other types of boundary conditions the factor be may 
depend on L. If one selects waves with other frequencies, which are peaked around 
instead of z/q = 3.3 mHz for the 5 minute oscillations, then the value of cr^ at i>i can be easily 
calculated from af « 1000 

With the time dependence the Doppler signal is 

dnimiR, t) = d^i^{R) exp {i{uJnit - anlm)} (36) 

The cross-correlation ipj"^ {Ri, R2, r) can be computed for the phase speed filtered Doppler 
signal d^j:^{R,9,(f),t) = J2nim-^pi-^^^ni)dnimiRy'^) ^ similar manner as for equation (fT2l) 
and involves the product of the projected line of sight Doppler signals at the two locations 
A (6*1,01) and B (6*2,^2) on the solar surface (See Figures E]and[HD. Here we have replaced 
animini{.R) ^"7 the Gaussian frequency function Giiuj), which models the amplitude of the 
solar modes. 

T 

ijl{Rr,R2,T) = iy dliR^,t)dliR2,t + T)dt (37) 



In evaluating this cross-correlation we need the summation over m of the spherical harmonics 
and its derivatives, which are given in Appendix C. 

Collecting all the sums over m we define the total sum (see Appendix C equations 
dlARi.R^) =Y.dliMdl,M = /oPKcosA) + f^^n(^ + /2^J^^ (38) 

m 

Substituting the expressions for P/(cos A) and its derivatives from Appendix B equations 
flFr|) and (jSSD, we obtain for large I, small A, such that LA is large 



dURi, ^2) = lV ^ [/o cos (la - ^) - If, sin (la - ^) - I'f, cos (lA - ^)] (39) 
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where L = + 1) and / = (L — |). Equation (|39|) can be written as 



Aj cos I sin I LA 



Ca/J^cos(LA-^ + C 

(40) 



ttLA 

where, the phase factor C(^i) 0i) ^2, ^2) = tsji~^ (Bi/Ai), Ai = (/o — 1^/2), Bi = lf\ and 
Ci = ^/ Af + Bf. The value of L outside the trigonometric functions is evaluated at Uo and 
is L{uo) = y- = The phase factor C, in general is a function of frequency Unu but since 
the solar oscillations are band-limited and their power is peaked at the frequency tUo, C is 
evaluated at the central frequency uJq- This also simplifies the evaluation of the integral, and 
a closed form expression can be derived. The phase shift C, is measured in radians but can 
be converted to units of time by dividing by ojo- 

The cross-correlation is 

i;jl{R,,R2,T) = Y,FpiL,U;r,l)GUcOnl)cOs{cUnlT)d:^l{Rl,R2) (41) 

n,l 

Substituting the expression for d^i{Ri, R2) from equation fHU]) into equation fHT]) and trans- 
forming the product of the cosines to a sum we obtain 



Vp LOnl 



hiuJnlT,C) + f-hii^nirX)] 



(42) 

where, f+h{uJnir, C) = cos {ujnir - LA + ^ - () and /_/,(u;„ir, C) = cos {u^it + LA - f + C) 
correspond to positive and negative time lags respectively. Evaluating the inner sum in 
equation fH2l) in a similar manner as in equation fl2T|) one obtains 



ipf^ (i?i, i?2, r, Vp) = Af^{6uj, 6f, r, Tg, Tp) cos |^/^(r - r/^J + -| (43) 
The shifted phase travel time r/^^ due to the horizontal component is therefore, 

rfpn =Vp + — (44) 

Summing equation (H3l) over phase velocities we get the final cross-correlation. 

^%{R^,R2,t) = Y,^^^fp{RuR2,r,Vp) (45) 
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In deriving this equation we used the asymptotic formula for the Legendre functions, which 
allowed us to obtain the explicit expression for the phase shift. We verified by direct nu- 
merical calculations of the cross-covariance functions that this approximation is sufficiently 
accurate. A comparison of the analytical formula in equation fllUl) for d^i{Ri,R2) with the 
direct numerical mode summation is shown in Figure [H This approximation holds when the 
travel times are much greater than the oscillation period of five minutes, which is the case 
for typical time-distance measurements. 




-180 -120 -60 60 120 180 



Fig. 1. — Comparison of function d^i{Ri, R2), calculated numerically by direct mode sum- 
mation given by the middle term in equation [38] (solid curve) and the analytical expression 
equation Ho] (dotted curve). In this example, the central ^1 point is located at 9i = 60° and 
01 = 30°, point R2 is located at a distance A = 8.4° from and the angle, 71 between the 
directions from Ri to the North pole and point R2 varies from —180° to 180°, as illustrated 
in Figure El Appendix D. 

From the equations (139|) and (HTj) we observe that the cross-correlation of the Doppler 
line of sight velocity involving the horizontal component of the displacement is a weighted 
sum of the phase speed filtered Gabor wavelets, with the weights depending on the location 
(6*1, 01, 02, 02) of the two points being cross-correlated, and the angular degree / of the modes. 
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This weighted sum can be conveniently combined incorporating a phase factor into the cosine 
term. 

From equation flH|) we see a shift in phase travel times due to the horizontal component 
in addition to that introduced by the phase speed filtering. In Figure [2], we see a variation 
in the travel time shift due to the horizontal component as we move point B (see Appendix 
E, Figure [H]) around the annulus by changing the angle 71 (see Appendix D, Figure [6]), cross 
correlating with the fixed point A. Figure E] illustrates the phase shift for the projection 
on the solar disk. The greatest phase shift changes happen for the outmost points in the 
direction from the disk center. 

Averaging the travel time shift over the annular angle 71, for different horizontal dis- 
tances A and distances from disk center rj, to simulate the observed mean travel times, we 
show in Figure H] that this time shift increases (becomes more negative) as we go away from 
disk center and also for increasing horizontal distances. The shifts are quite appreci able, on 



the o rder of a few seconds, and have been observed in the real data from SOHO/MDI (IDuvall 



20031 ) . At disk center C, (see Appendix E, Figure [7]) which according to our convention is at 
{y, z) = (0, 0), which is the origin of the YZ plane perpendicular to the line of sight direction 
Cx- In polar coordinates the disk center is specified by {6, (p) = (^, 0), at which the horizontal 
component is zero in equation flMl) . and the radial component is maximum. We see from 
Figure [21 when point A is at disk center, the time shift is independent of 71, moreover there 
will always be a small travel time shift for this case since point B is at a distance A from 
the disk center, and it has a non-zero horizontal component that when cross-correlated with 
the radial component of the signal at point A, which is at disk center, leads to a travel time 
shift. 



Dispersive and Non-Dispersive Wave packets: Relation to the Gabor 
Wavelet without Phase Speed Filtering 



A dispe rsive medium i s where the pha se velocity Vp is quite different from the group 
velocity Ug ( lLighthilllll978l : IWhithamlll974l ). Due to this for a fixed distance A, the group 
and phase travel times are different. In the solar case the wave packets that travel short 
distances, probe the outer layers of the Sun, which are quite dispersive due to the presence 
of large gradients. Hence the wave packets got by the cross-correlation in equations fl2^ and 
( 127]) are distorted as they travel through the medium. Moreover, the approximations made 
by retaining only the first order terms in the Taylor expansion of the dispersion relation in 
equation ([8]) may be inadequate to completely model the physical effects associated with 
small distances, hence care must be exercised in interpreting the results got by using these 
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Fig. 2. — Phase shift (Y axis) ^ in seconds as a function of angle 71 (X axis) in degrees. 
The top most row corresponds to (pi = —10 degrees, the next rows are in increments of 
10 degrees. The left most column corresponds to 61 = 90 degrees, the other columns are 
in decrements of 10 degrees. With this convention the left most plot of the second row 
corresponds to {61, (pi) = (90°, 0°), for the annulus central point A, which is at disk center 
C. The star (*), plus (+) and cross (x) lines respectively in each plot correspond to A = 
9.84, 15.36 and 20.64 degrees, and Tp = 58.89, 69.24 and 77.01 minutes respectively. In this 
figure there is no shift due to phase speed filtering. 



wavelets for small distances. 

On the other hand, wave packets that travel large angular distances A, probe the deeper 
layers of the Sun, that are weakly dispersive, and hence undergo less distortion. It is quite 
interesting to note that the Gabor wavelet in equation (13T]) . in the absence of the phase speed 
filter is got in the non- dispersive limit of equations (126!) and ( !27|) when Rg = Rp = 0, that is 
Tg = Tp = Tph, for the intensity or line depth signal which has only the radial component. This 
is quite intuitive since in this limit there is no phase speed filtering, and Fp{L,u) = 1, hence 
the phase speed filtered signal in equation (ITTi) reduces to the frequency band-limited signal 
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X 



Fig. 3. — Illustration of the phase shift variations for a projection on the solar disk. The 
travel times are measured between the central point, located at 9i = 60 degrees and 0i = 30 
degrees, and surrounding points at distance A = 8.4 degrees and corresponding Tp = 55.59 
minutes. The size of the symbols is proportional to the magnitude of the phase shifts which 
varies from -0.07 to -0.12 radians (corresponding to travel time shifts from -3.5 seconds to 
-5.5 seconds). In this Figure there is no phase speed filtering. 



in equation This is consisted both physically and mathematically. To form wave packets 
the medium need not be dispersive. Waves in a range of frequencies and wave lengths can 
interfere to form a wave packet. The wave packets that travel through a dispersive medium 
are distorted, where as those through a non- dispersive medium are not distorted, but retain 
their original shape. 

The dimensionless quantities Rp and Rg measure the amount of deviation from non- 
dispersiveness. Hence to fit the observed temporal cross-correlation for the intensity or line 
depth signal where, Rg and Rp are small, there will be little error in using the Gabor wavelet 
of equation fl^ . since t he new Gabor wavelet of equations (12^ and (1271) is very close to 
the old Gabor wavelet ( Kosovichev fc Duvall 1997 : Giles 19991 ) in this non-dispersive or 
weakly dispersive limit for the line depth or scalar intensity signal that have only the radial 
component. 
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Fig. 4. — The mean phase shift in seconds (Y axis) got by averaging ^ over the annular 
angle 71 as a function of ratio with respect to solar radius of distance of the central point 
A from disc center C, ^ (X axis). The star (*), plus (+) and cross (x) lines respectively 
correspond to A = 9.84, 15.36 and 20.64 degrees, and Tp = 58.89, 69.24 and 77.01 minutes 
respectively. We observe a maximum travel time shift of about 6 seconds. This plot has 
no phase speed filtering, whose effect can be removed by choosing Vph appropriately (see 
Sec.2.1). 



4. Conclusion 



In this paper wc introduce a mathematical formalism to deal with the horizontal com- 
ponent of the displacement, which is used to derive a new formula to estimate travel times 
in time-distance helioseismology. It models the observed cross-correlation, by considering 
both the radial and horizontal components of the displacement, and also includes the phase 
speed filter on the frequency band-limited solar oscillations. The form of the Gabor wavelet 
is retained by the present formula, and in the non-dispersive limit for the line depth or inten- 
sity observations it reduces to the previously used Gabor wavelet that was derived without a 
phase speed filter, considering only the radial component. Including the horizontal compo- 
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nent, the cross-correlation is a weighted sum of phase speed filtered Gabor wavelets, with the 
weights depending on the angular position of the two points being cross-correlated, and the 
angular degree of the modes. Due to the phase speed filtering and inclusion of the horizontal 
component of the displacement, shifts in phase travel times are observed. In particular, 
the horizontal component induces a systematic shift in travel times, whose absolute value 
increases as we go away from disk center and also for large horizontal distances. This will 
in turn have implications in the inversions to study flows and the sub-surface properties in 
the Sun and also to formulate other local helioseismic approaches that involve the horizontal 
component. 
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6. Appendix A 

We define the spherical coordinate system used in this paper (see Figure [5]). 

Consider fixed cartesian unit vectors {e^, Cy, e^) and the corresponding spherical unit 
vectors {ir, ee, e^) that vary with the angles {9, 0). The equations relating these sets of unit 
vectors are 

= Cr sin 9 cos 4> + eg cos 9 cos — e<^ sin (46) 
iy = Cr sin 6^ sin + ee cos 6^ sin + cos (47) 
e-z = Cr COS 9 — eg sin 9 (48) 

Any point (x, y, z) in spherical coordinates (r, ^, 0) is x = r sin 9 cos (p, y = r sin 9 sin 
and z = rcos9. Consider two points A (^i, 0i) and B {92 , 02) on the solar surface r = R 
(see Figure El Appendix D), with coordinates -Ri = {xi,yi,Zi) and R2 = (2:2,1/2,-22), the 



This preprint was prepared with the AAS M^TjX macros v5.2. 
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Fig. 5. — Spherical Coordinate System. The hne of sight direction is along e^. 

angular distance A between them can be calculated by taking the dot product of the unit 
vectors at these two locations and is given by 

_ (-Ri • R2) _ X1X2 + yiy2 + Z1Z2 , . 

cos A- ^2 - R2 l^yj 

where, {xi,yi,zi) = (i? sin 6*1 cos 0i, i? sin 6^1 sin 0i, i? cos 6*1) and 

{x2,y2,Z2) = (-R sin ^2 COS 02, -R sin 6*2 sin 02, -R COS 6*2) (see Figure [6l Appendix D). Substitut- 
ing this into equation (H9l) we obtain 



cos A = cos 9i cos 62 + sin 9i sin 62 cos(0i — 02) (50) 



7. Appendix B 

We derive various approximations f or P/(cos A) an d its derivatives for large I, small A 
and large LA following the approach of (|Jacksonlll999l ). 



Pi(cosA) ^ Jo 



{21 + 1) sin 



A 



ttLA 



cos ( LA 



(51) 
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From the recursion relation 



g-^^ ^-P^cosA) _^cosAPKcosA) + ^^^-i(cosA) =0 (52) 
dA 

For small A equation (1521) becomes, 

^dPi{^A)_ _ ^p^^^^g^^ ^ /P;_i(cos A) = (53) 

Differentiating equation flS^ with respect to A gives 

d'^PAcosA) , dPi(cosA) dPi^i(cosA) 

sin A -i— — - - / - 1 cos A — ^- ^ + / — ^ + / sin AP, cos A = 54 

dA'' dA dA 

For small A and large /, equation fl54|) becomes, 

^d^P,(oos A) _ ^dP.(cosA) ^ ^dP,_.(cos A) ^ ^ ^ ^^^^ 

(JjL^ UjL^ (JjL^ 

Now for small A, large / and large LA the various approximations are after using the recursion 
relations from equations (153|) and (l55ll 



P^cosA) ^ ^/-|^cos(lA-^) (56) 



dPi{cosA) 



dA V ttLA 



sin (la - ^) (57) 



d^Pi(cosA) ,o / 2 / . 7r\ 



c/A2 V ttLA V 4 



8. Appendix C 



Using the addition theorem of the derivatives of spherical harmonics (jWinch and Roberts 



19951 ). we derive the expressions for the terms in equation (l38l) . 
The angular distance A is, 

cos A = cos 9i cos 02 + sin 9i sin 02 cos(0i — 02) (59) 



Differentiating equation fl59l) with respect to 0i, 02, 0i and (j)2 we obtain 
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sin 6i cos 62 — cos 61 sin 62 cos(0i — 
96'i sin A 



COS71 



(60) 



dA cos 9i sin 6^2 — sin 9i cos 6*2 cos(0i — 02 ) 



sin A 

dA sin ^1 sin 62 sm{(f)i — 02) 
c}0i sin A 

dA sin 6'i sin 62 sin(0i — 02) 
(909 sin A 



cos 72 



sin 61 sin 71 



sin 62 sin 72 



(61) 



(62) 



(63) 



The e xpressions in equations ([60l) - (!63ll can be equated to sine and cosine functions of 
7i and 72 ( IWinch and Roberts! llQQSh. which a re shown in Figure El Appendix D using the 
results of spherical trigonometry ( ISmartI 119771 ) . 



According to the addition theorem of spherical harmonics (lJacksodll999l ) 

Snl = 5^rr(^l,0l)Fr(^2,02) = ^^^P^COSA) (64) 



For notational convenience we denote Ni 
respect to 62 and applying the chain rule leads to 



(21+1) 



. Differentiating equation (16^ with 



72, V2) 



'li VI) 



89, 



N, 



dPi{cosA) OA 
dA 



(65) 



Differentiating equation flMl) with respect to 02 and applying the chain rule leads to 



1, <P\) TT, 



dPi{cosA) dA 
dA 



(66) 



Differentiating equation fl64l) with respect to 61 and applying the chain rule leads to 



Or 



dA d9i 



(67) 



Differentiating equation (167]) with respect to 62 and applying the product and chain rule 



gives 



dYr{ei,<Pi)dYl 



72, (P2 



002 



N, 



dAdAd^Pi{cosA) ^ d f dA\ dPi{cos A) 



ddi 802 dA2 



062 ydOi 



dA 



(68) 
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Differentiating equation fl67l) with respect to 02 and applying the product and chain 
rule gives 



^ dYr{eu<Pi)dYr{e2A2) _ 

^ dOi d(t)2 ' 



rf2Pi(cosA) ^ d ^ dA\ dPi{cos A) 



d9i 



2 ydej dA 



(69) 



Differentiating equation (1641) with respect to (pi and applying the chain rule leads to 



<j)r 



h,<l>l)ym.n , , ,, rfPKcOsA) dA 



(70) 



Differentiating equation fjTOj) with respect to 02 and applying the product and chain rule 



gives 



_^ dYrieu(Pi)dYr{e2A2) _ ,r 

'nl - 2^ JU. - 



OAdAd^PiicosA) d f dA \ dPi{cos A) 



dA 



(71) 



Differentiating equation (ITOl) with respect to 02 and applying the product and chain 
rule gives 



dYr{0i,<Pi)dYr{e2,<P2) 



N, 



OA OA d^Pi (cos A) d fdA \ dPi (cos A) 



h d(j)2 rfA2 



dA 



(72) 



Now using equations ([31]), flMl) - fl721) we obtain for the various terms in equation 



4z = sin 6*1 cos 01 sin 6^2 cos02<'? = /[q (6*1, 0i, 6*2, 02)^«(cos A) 



(73) 



For notational convenience we denote the angular coordinates of points Ri and R2 by 

ai2 = (6'l,01,6'2,02) 



sll = Pnl{R) sin 6^1 cos 01 cos 6*2 COS(f)2Sn^i = f2l{^12 



dPi{cosA) 
' dA 



(74) 



3 p ■ a A fr<Pf^ dPijcosA) 
Sni = -l3ni[R) sm^i COS 01— J = /3r(ai2) — 

er _ fen-t .dPi{cosA) 



Snl = l3nl{R) COS 6*1 COS 0i siu 6*2 COS 02S J = /4[(ai2 



dA 



(75) 
(76) 
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Sni = -[PniiR)] COS 01 cos0i^t^sJ = fQ^{au) + fdiau) (7J 



7 Q fjD\^^^^^ ■ a A 'f'r f^ri^ dPi{cosA) 
Snl = -Pni{R)^^sme2Cos(f)2sli = fji [au) (79) 



8 in /DM2 sin 01 n<t>9f^ d'^PijcosA) dPijcosA) 

Snl = -[MR)] COS ^2 COS02<, = f^2 M + ^Sl M ^ (80) 



9 /mi2Wisin02 <^^_ d^PijcosA) dPijcos A) 

Where we evaluate the derivatives of A with respect to 9^, 6 2, 0i and 02 using equa- 



tions flBn]) - fl^ and the equations from ( Winch and RobertillQQsl ). we obtain the following 
expressions 



/lo (^12) = Ni sin 6*1 cos 0i sin O2 cos 02 
/21 (^12) = Nil3ni{R) sin 6^1 cos0i cos 6^2 cos 02 cos 72 

/af (012) = Ni(3ni{R) sin di cos 0i sin 02 sin 72 
/41 (012) = Ni(3ni{R) cos 9i cos 01 sin 62 cos 02 cos 71 

ftlidu) = -A^z[/?n«(-R)]^COs6'iCOS0iCOs6'2COS02 



sin A 

/52(^12) = Ni[l3nl{R)]'^ COS 61 COS 0i COS 6*2 COS 02 COS 7i COS 72 

ftiidu) = Ni[(3ni{R)]'^ COS 611 COS 01 sin 02- 



sin A 

/62^(ai2) = A^Z [/3n« (R)? COS 9i COS 0i slu 02 COS 7i slu 72 

ff[ = -Ni(3ni{R) sin 01 sin 62 cos 02 sin 71 



82) 
83) 
84) 
85) 

86) 

87) 

88) 

89) 
90) 
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fi?{Su) = -Ni [(3ni (i?)] ^in 01 cos 9^ cos 02 ^^^^ 

fwi^i^) = -Ni[f3ni{R)]^ sin (pi cos 02 cos (j)2 sin 'y I cos '^2 (92) 

/'/i'/'^- ^ Arr/? ^dm2 • • cos 71 cos 72 , . 

/91 (012) = A^z Pni(^) sm0ism02 — (93) 

sm A 

/92^(«i2) = -Ni[Pni{R)]^ sin 01 sin 02 sin7i sin 72 (94) 

In computing the cross correlation in equation (HTl) . we need to evaluate d^i{Ri,R2) in 
equation (l38ll as 

dURuR2) = J2<ln.{Rl)dnlrniR2) = E = /o^^COS A) + ^ /-P^^OS A) ^ rf^P^^OS A) 

(95) 

where, 

/o(ai2) = /[o' (96) 

/i(ai2) = /2t + fst + fll + /5T + fit + /rt + ft! + /s? (97) 

f2{d,2) = ft^ + flt + ft2+ftt (98) 



9. Appendix D 

In this appendix we derive the angular coordinates (6^2,02) of point B on the annulus, 
which is at a fixed annular angle 71 and a given angular distance A from the fixed point A 
with angular coordinates (6*1,01) at the center of the annulus. We use equations (!59|) and 
f l60p from Appendix C to express the angular coordinates of B in terms of the fixed angles 
(^1, 01, 7i, A). The geometry for the calculation is shown in Figure 6. 

cos A = cos 61 cos 62 + sin 61 sin 62 cos(0i — 02) (99) 
sin A cos 7i = sin 6*1 cos 62 — cos 61 sin 62 cos(0i — 02) (100) 

We now consider the case when the two points are close together, which is true for small 
distances A. They are related as 6*2 = ^1 + 56i and 02 = 0i + ^0i, for small increments 56*1 
and 501. The equations (p9|) and fllOOp can be expanded in terms of the increments, and 
using the standard trigonometric approximations for small angles sin 9^9 and cos ^ ~ 1 — y 
we obtain. 
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Fig. 6. — The Spherical Triangle showing A, which is the angular separation of the great 
circle arc between points A and B on the solar surface. Here, N is the North pole. 



861 = A 



1 - 



sin^ 7i sin^ 9i 



sin ^1 — "2 cos 7i sin 26i 
A sin 7i 
[sin^ 9i — ^ cos 7i sin 29i^ ^ 



(101) 
(102) 



For fixed A and 71 we use equations f llOip and f ll02p to calculate (6*2, ^2) for small 
angular separation. For larger separations we can solve the linear equations fl^I?]) and fllOOl) 
to calculate 62, and equation (l62l) to compute 02, and equations (16T1) and (163|) to calculate 
cos 72 and sin 72 respectively. 



10. Appendix E 



In this appendix we consider the orthographic projection shown in Figure [TJ where the 
points on the surface of the Sun are projected onto the YZ plane P, which is perpendicular 
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Fig. 7. — Orthographic projection. The point C is the disk center and rj measures the 
distance from C to any projected point (Y,Z) on the orthographic plane P, which is at 
X — R. Any point A' on the solar surface is specified by (^i, 0i) or the angle 0i and the 
distance 77. 

to the line of sight direction e^. The origin O is at the center of the Sun of radius R. We 
denote a point A' on the solar surface as {x, y, z) with axis xyz and as (F, Z) the coordinate 
of the same point projected as A on the projection plane P with axis YZ . The x axis is 
perpendicular the the projection plane. The projection is oriented with north along +y axis 
and east along the +z axis. The projection is centered at co-latitude of 90° and longitude 
of 0°. With this convention, the north pole N has coordinates (0,i?, 0), the south pole S is 
(0, —R, 0), and the point of tangency C between the sphere and the plane is (i?, 0, 0). The 
point C is also called the disk center. The projection plane \s x — R. The coordinates of 
any point A' with co-latitude 9i and longitude 0i on the solar surface are 

{x,y,z) = (i? sin 6*1 cos 01, sin ^1 sin 01, cos 6*1) (103) 

The equation of the sphere is x'^ + y'^ + z'^ — R^. The coordinates of the same point on the 
plane are (Y, Z) — {y,z). The point of tangency C between the sphere and the plane is the 
disk center and has coordinates (Y, Z) — (0, 0) on the plane. Hence, the distance from the 




Fig. 8. — Projected Disc. This shows the annulus D/^ of radius A around the point A (^i, (pi) 
on the solar surface. The cross correlation is computed between the annulus center A and a 
point B {621 (P2) on the solar surface, that has an angular separation of A from A. The point 
A is at a distance r] from the disk center C, {9^, (pc) = (f , 0). 

disk center C (0, 0) to any point A [Y, Z) on the plane P is denoted by rj, which is fixed 
for each computation of the cross correlation. Any point A' on the solar surface and hence 
the corresponding projection A on the plane P can be specified by (6'i,0i) or alternatively 
by the two parameters t] and the angle 0i. In Figure El the center to annulus geometry is 
shown with the fixed center point A {61, (pi), which is at a distance 77 from disk center C 
(y, Z) = (0, 0), and a point B (6*2, ^2) is on the annulus of radius A. 

7] = {Y^ + Z^)^ = {y^ + z")^ = {R^ - = R{1 - sin^ 61 cos^ (104) 

The angular coordinates of the disk center C, are {dc,(l)c) = (f;0). When the point A 
is at the disk center (6'i,0i) = (90°, 0°), and hence the distance r] = 0. For fixed values of 
(A,7i), together with the coordinates (6^1, 0i) of the annulus center A, we can compute the 
coordinates {62, 02 ) of point B on the annulus by solving the linear system of equations fl99l) 
and fllOOl) to get 62 and solving equation fl62l) to compute 02- 



